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Abstract 

Using the Bose-Einstein correlations (BEC) implemented in PYTHIA we investigated a possibility of the 
CDF experiment at the Tevatron to see the two-particle correlations in the final state of interactions. 
The approach based on quantum field theory at finite temperature was applied to the ALEPH data at 
LEP, and the BEC important parameters were retrieved. 



1 Introduction 

Over the past few decades, a considerable number of studies have been done on the phenomena of 
multi-particle correlations induced by collisions between elementary particles. It is understood that the 
studies of correlations between produced particles, the effects of coherence and chaoticity, an estimation 
of particle emitting source size play an important role in high energy physics. 

In papers [1,2], there were established and developed the master equations in the form of the field operator 
evolution equation (Langevin-like [3]), which allows one to gain a better understanding of the structure 
of the emitting source. 

The shapes of Bose-Einstein correlation (BEC) function were established in the LEP experiments ALEPH 
[4], DELPHI [5] and OPAL [6], and ZEUS Collaboration at HERA [7], which also indicated a dependence 
of the measured correlation radius on the hadron (tt, K) mass. 

There is still no definite explanation of the above mentioned mass dependence. The aim of this paper 
is, firstly, to carry out the numerical analysis applied to the LEP data (e.g., ALEPH Collaboration) on 
two-particle BEC in the framework of quantum field theory at finite temperature [QFTp) approach [1,2], 
and secondly, to give the proposal on BEC studies at the Tevatron energies. 

We simulated pp annihilations at the Tevatron energies (1.96 TeV). The following relation was used to 
retrieve the BEC function C2{Q) from simulated data: 

C,{Q) = N^^{Q)lN^^f{Q), (1) 



where N^^{Q) is a number of like-sign pairs with four momentum difference Q = y (-Pi ^ Pa)^, [i = 

1, 2) are four-momenta of produced particles. N^'^^ is the number of pairs without BEC. The reference 
sample can be created in two different ways. One can use the unlike-sign pairs from the same event or 
take the like-sign pairs from different events. Both approaches were employed in the present work. 
The correlation procedure in PYTHIA is based on the following assumptions. First of all, the Qij = 
\J {Pi + Pi)^ ^ value is evaluated for each pairs of produced particles with the mass m. Then the 
shifted Qi,- is found as a solution of the following equation 





using the chosen parameterization of C2{Q) function. 



(2) 



2 Experiment ALEPH vs. QFT^ model of BEC. 

In general, the shape of BEC C2 [Q) function is model dependent. The most simple form of Goldhaber-like 
parametrization for C2((3) [8] 



C2(Q)-eW(l + Ae 



)(i + eQ), 



(3) 



has been used for data fitting in LEP and HERA experiments. In formula (3), £,{N) is the normalization 
factor, depending on multiplicity N , A is so-called the coherence strength factor, meaning A = 1 for 
fully incoherent and A = for fully coherent sources; the parameter R is interpreted as a radius of the 
particle source, assumed to be spherical in this parameterisation [8,9]. The linear term in (3) accounts 
for long-range correlations outside the region of BEC. We use the C2{Q) function within QFTp approach 
[1,2] in the form: 

C2{Q)=aN) l + -^^^) + -^^{q) , (4) 
where the multiphcity N depending factor S,{N) is equal to: ^{N) ^ {N {N ~ 1)) / {Nf. 
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Figure 1: The one-dimensional correlation function constructed using the ALEPH data with the mixed 
reference sample (left), C™"(Q), and with the unhke-sign reference sample (right), 6*2"'*'^'^ (Q). The data 
are fitted with Eqs. (3), (4) and (5). 







C+-(Q) 






0.948 ± 0.001 


0.936 ± 0.001 




X 


0.362 ± 0.006 


0.438 ± 0.006 




R{fm) 


0.528 ± 0.005 


0.777 ± 0.007 


(3) 


e 


(0.768 ± 0.01)10-2 


(0.905 ± 0.002)10-2 




xVndf 


513 / 94 


432 / 72 




m) 


0.993 ± 0.001 


0.9871 ± 0.0007 




a 


4.32 ± 0.08 


3.28 ± 0.08 




n{fm) 


0.827 ± 0.008 


1.21 ± 0.01 


(4) 


xVndf 


834.5 / 95 


1108 / 75 




m) 


0.946 ± 0.003 


0.938 ± 0.002 




a 


3.85 ± 0.07 


2.92 ± 0.06 




R{fm) 


0.713 ± 0.009 


1.04 ± 0.01 


(5) 


rf(fm) 


(0.79 ± 0.05)10-2 


(0.87 ± 0.04)10-2 




xVndf 


483.7 / 94 


465.8 / 74 





Tabic 1 : Comparison of the fit output parameters obtained with different C2 (Q) functions. The ALEPH 
data were fitted with Eqs. (3), (4) and (5). 



The important parameter a summarizes our knowledge of other than space-time characteristics of 
the particle emitting source. Q{q) is a Fourier image of hadronizing source given by ^{q) = e~'^ 'ji'n) = 



7(n), where j{n) 



n{uj)n(uj') ' 



and 



n(a;) = n(a;, /3) 



g(a>-M)/3 _ 1 ' 



where n{oj, (3) is the number of particles with the energy co in the thermal bath with temperature T = ^. 
For a more complete analysis we use the modified formula for C2(Q): 



C2(Q)=^(iV) 



2a 



(1 + a)' 



1 



(1 + r/Q), 



(5) 



where the consequence of the Bogolyubov's principle of weakening of correlations at large distances was 
not ignored; r/ is the parameter of weakening of correlations. 

The correlation functions, C2"'*'^'^(Q) and C™"((5), obtained using the ALEPH data, are shown in Fig. 
1. In Table 1., we present the main retrieved parameters containing in formulae (3), (4), and (5) using 
the ALEPH data. 



3 PYTHIA Simulations 

We fit the simulated data from PYTHIA with the most common C2{Q) function (3) at e = 0. In the 

first step, we reproduced the parameters A and R, that had been put into PYTHIA at the beginning of 
simulation. We have found out there are non-vanishing correlations between 7r° mesons even if we turn 
off the BEC. Therefore in the next analysis we take into account only tt"*" and 7t~ mesons. 

Fig. 2 demonstrates how the output parameter R behaves as func:tion of input parameters. The 
green area around = 1.4/m (right plots) is the area of the best reproduction of the parameters. It 
means, that the difference between input and output value of parameter AR = Rout — Rin in this area is 
approximately zero. The red area in right bottom side corresponds to Ai? = 0.2 fm, and the blue area in 
top left side to AR = — 0.5/m. The output parameter A has no a good reproducing area. Each output 
value of the parameter A was smaller than the input one in the range 0.2 up to 0.5. 




Figure 2: The dependence of the output parameter Rout on the input parameters i?i„ and A (left). The 
different colors correspond to different values of Ai? = Rout — Rin (right). 

4 Two-particle BEC for CDF 

In the next step, we have investigated a possibility to study the two-particle correlations in the CDF 
experiment. We blurred all components of momenta of each particle in accord with the CDF detector 
resolutions [10]. The PYTHIA simulations show us, that a lot of pions have energy too small for reaching 
the hadron calorimeter. Therefore we get new particle energies using blurred momenta and particles 
masses (Fig. 3), employing, for the blurring procedure, the realistic CDF tracking system resolution. 
It is evident that taking into account realistic CDF detector resolutions we are able to observe the BEC 
effects, although the C2{Q) function is less profound in this case. 

In Fig 4. one can see a decrease of the output parameter Rout, taken as a function of the input parameter 
Rim with increasing cut on energy. For suppression of the decrease it is therefore important to register 
also low energy particles. 



5 Summary 

Based on QFTp we carried out the numerical analysis of the ALEPH data applied to two-particle BEC. 
As a result, the important parameters of C2{Q) function {R,a,£,{N),rf) are retrieved. We investigated 
the possibility of the CDF detector at the Tevatron to explore the two-particle correlations in the final 
state. As a result, there will be possible to observe the BEC effects in case of charged pions although 
the 6*2 (Q) function may be slightly suppressed. The importance for measuring of low-energy particles 
for BEC studies at the CDF is established. 
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Figure 3: The dots (black curve) correspond to 
the C2{Q) function obtained from PYTHIA with no 
blurring and the squares (red curve) correspond to 
the C2{Q) with blurring. 



Figure 4: The dependence of output parameter R on 
input parameter R and different energy cuts; A is set 
to 0.4. 
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